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<N ' Abstract 

The paper deals with a semilinear integrodiffcrential equation that character- 
izes several dissipative models of Viscoelasticity, Biology and Superconductivity. 
Cn| . The initial - boundary problem with Neumann conditions is analyzed. When 

the source term F is a linear function, then the explicit solution is obtained. 
When F is non linear, some results on existence, uniqueness and a priori es- 
timates are deduced. As example of physical model the reaction - diffusion 
system of Fitzhugh Nagumo is considered. 
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^ ■ 1 Introduction 



Let consider a function u = u(x,t), where x is a direction of propagation and t is 



£Nj ■ the time, and let 

.>: t 

(1.1) Cu = u t -eu xx + au + b e -/3( '" r) u(x,t) dr = F(x, t, u(x,t)) 

j_j \ Jo 

with a, b, e, /3 positive constants. 

The equation (1.1) describes the evolution of several physical models as motions 
of viscoelastic fluids or solids [3,6,17]; heat conduction at low temperature [13], 
sound propagation in viscous gases [12]. Other two specifical examples for the integro 
differential equation (1.1) are related to biological models and superconductivity. 
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As for the biological phenomena, a well known reaction diffusion model is given 
by the Fitzhugh - Nagumo system (FHN) [9,14,19]: 



du d 2 u 



dt dx 

;i.2) 

— = bu - pv. 
dt 



In this case the function f(u) is: 

(1.3) f(u) = — au + (f(u) with ip = v 2 (a + 1 — u) (0 < a < 1) 

As for the variables u and v , u(x, t) represents a membrane potential of a 
nerve axon at distance x and time t, and v(x,t) is a recovery variable that models 
the transmembrane current. 

If vq represents the initial value of v, the system (1.2) can be given the form 
(1.1) with 

(1.4) F(x,t,u) = <p(u) - v (x)e~ pt . 

Moreover, equation (1.1) occurs also in superconductivity to describe the Joseph- 
son tunnel effects in junctions. In this case the unknown u denotes the difference 
between the phases of the wave functions of the two superconductors and the dif- 
ferential equation is: 

(1.5) eu xxt - u tt + u xx - au t = sinu +7 

where 7 is a constant forcing term that is proportional to a bias current. The 
e — term and the a — term account for the dissipative normal electron current flow 
along and across the junction, respectively [1, 18]. 

From (1.1) one obtains the equation (1.5) as soon as one assumes 

(1.6) a = a-- b=-- p =- 

e e e 

and F is such that 



(1.7) F(x,t,u) = - e- = ( *- r) [7 + senu(x,T)} dr. 

Jo 

As (1.6) show, in the superconductive case the constants a, b could be negative 
too. 

The explicit fundamental solution Kq(x, t) of the operator C defined in (1.1) has 
been already determined in [7] together with numerous basic properties. When F is 
a linear function, by means of Ko(x,t) it is possible to obtain the explicit solution 
of both the Neumann and the Dirichelet problem for (1.1). When F is non linear, 
an appropriate analysis of the integro differential equation implies results on the 
existence, uniqueness and a priori estimates of the solution. These results will be 
applied to (FHN) system. 

2 Statement of the problem and tranform solution 

If T is an arbitrary positive constant and 

Q T = { (x, t) : < x < L ; < t < T }, 

let (-P/v) the following Neumann initial - boundary value problem related to equation 
(1.1): 



(2.1) 



ut — eu xx + au + b f Q e &" T > u(x,t) dr = F(x,t,u) (x,t)E&T 
-u(x,0) = uo(x) x £ [0, L], 

I u x (0,t) = Mt) u x {L,t) = Mt) 0<t<T. 



In excitable systems this problem occurs when two-species reaction diffusion system 
is subjected to flux boundary condition [15]. The same conditions are present in 
case of pacemakers [11]. Neumann conditions are applied also to study distributed 
FHN system [16]. 

In superconductivity, instead, (-Pat) problem can be referred to the boundary speci- 
fication of the magnetic field [2,8, 10]. 

When F = f(x, t) is a linear function, the problem (Pjv) can be solved by Laplace 
transform with respect to t. 



If 



(2.2) 



u(x,s) = / °° e st u(x,t)dt f(x,s) = / °° e st f (x,t)dt , 
k & ( s ) = io°° e ~" V>i (*) eft (^ = 1,2), 



one deduces the following transform (Pn) problem: 



Hirr-rr. (Jj 



(2.3) 



±[f( x ,s) + u (x)} 



u x (0,s) = i>i(s) u x (L,s) = 4> 2 (s), 



where a = s + a + -. Letting a = a Is, and considering the following 

s + p 
function 



(2.4) 



%{V,&) 



cosh [a (L — y)\ 
2 e a sinh (a L) 



1 



2 ^ a 



e v^ 



* + £ 



n=l 



2n.L-f-T/ 2nL — y 

e v^ + e ^ 



the formal solution u(x,s) of the problem (Pn) can be given the form: 



(2.5) u(x,s)= [0 (\ x -£\,s)+ d (\x + S\, s)][u U) + ftt,s)]dt 

Jo 

- 2 £ fa (s) § (X, S) + 2 £ 1p2 (s) 9 (L -x,s). 



3 Explicit solution in the linear case and asymptotic 
properties 



The fundamental solution Kq(x, t) of the linear operator C defined in (1.1) has been 
already obtained in [7] and it is: 



(3.1) K (r,t) 



2 
1 



2V^ 



e -Tt~ at n f t e ~*5~ aV 



y/t Jo Vt-y 



Vb , e-^ t ^yU 1 (2^by(t-y))dy 



where r = \x\ j \fe and J n {z) denotes the Bessel function of first kind. 

More, the following theorems have been proved in [7]: 

Teorema 3.1. In the half-plane Ke s > max( — a, — B ) the Laplace integral Ct Ko(r, t) 
converges absolutely for all r > 0, and it results: 



roc e~ ra 

(3.2) CtK = / e~ st K (r,t)dt = —^. 

Jo 2 Veer 



Teorema 3.2. The function Kq has the same basic properties of the fundamental 
solution of the heat equation, that is : 

i) K (x,t) eC°° for t>0, x G». 

ii) For fixed t > 0, Kq and its derivatives are vanishing esponentially fast as 
\x\ tends to infinity. 

Hi) For any fixed 5 > 0, uniformly for all \x\ > 5, it results: 



(3.3) lim K x,t) = 0, 



iv) For t > 0, it is C K = 0. 

Moreover, if u = min(a, (3) and one puts 



e"^ - e~ at 1 r- a + B 

(3.4) E(t) = 1 > 0, Bo = - + irVb 



a - 8 ' a 2(a/3) 3 / 2 ' 

then the following estimates hold [7]: 

(3.5) \K \ < -^= [e~ at + btE(t)}; f dr f \K (x - £,t)\ d£ < A, 

2\/TTet Jo Jsft 



(3.6) 



/ \K (x-£,t)\ d£ < e~ at + VbTTte' 



w t 



In order to obtain the inverse formulae for (2.5), let apply (3.2) to (2.4). Then 
one deduces the following function symilar to theta functions: 



(3.7) 9 (x,t) = K (x,t) + ]T [K (x +2nL, t) + K (x - 2nL, t)\ 

n=l 
oo 



As consequence, by (2.5), the explicit solution of the linear problem (-P/v) where 
F = f(x,t) is: 

(3.8) u(x,t)= [ [e o (\x-S\,t)+0 o (x + Z,t)]uo(Od£ + 

Jo 

-2s 6 (x, t-r) ^i(r) dr + 2 e / 9 (L-x,t-T) tjj 2 (t) dr 
Jo Jo 

+ [ dr [ [e (\x-$\,t-T) + 6 (x + Z, t-r)} f(£, T ) d£. 

Jo Jo 



Owing to the basic properties of Ko(x,t), it is easy to deduce the following 
theorem: 

Teorema 3.3. When the linear source f(x,t) is continuous in Sly and the initial 
boundary data Uq(x), ipi (t) (i = 1,2) are continuous, then problem (P/v) admits a 
unique regular solution u(x,t) given by (3.8). I 

As consequence of the properties of fundamental solution Ko(x,t), various esti- 
mates for u, ut, u x ... could be obtained. 

As an example, let evaluate the asymptotic properties of the terms caused by 
the initial datum uq(x) and the source f(x,t). If 

|| ito ||= sup \u (x)\, ll/H = sup | f(x, t) |, 

o<x<l n T 

it results: 



Teorema 3.4. When ipi = (i = 1,2), the solution (3.8) of (P/v), /or Zarge t, 
verifies the following estimate: 



(3.9) \u(x,t)\ < 2 [\\f\\ /3 + || i*o || (1 + V^Trt) e" 

where to = min (a, /3) and (3q is defined by (3.4)2- 
Proof: Properties of Ko(x,t) imply that: 



■u t 



(3.10) 



\ 6 (\x - £\, t) d* < V / \K (\x-t + 2nL\,t)\d£ 

JO „_*"L JO 



= V / |/ro(y,t)|dy < / |-Ko(y,t)|dy. 

n=-oo ^+(2n-l)£ 73? 

So, appling properties (4.10)2 and (3.6) to (3.8), the estimate (3.9) follows. I 

4 The Fitzhugh - Nagumo model. A priori estimates 

Consider now the non linear case of the (FHN) model defined by (1.2). By means of 
the previous results we are able to obtain integral equations for the two components 
(u, v) in terms of the data. All this implies the qualitative analysis of the solution 
together with a priori estimates. 

At first let us observe that by (1.2)2 one has: 

(4.1) v = voe~ pt +b [ e-^ t - T ^u(x,T)dT 

Jo 

and this formula, together with (1.4) require the presence of the following convolu- 
tions: 

(4.2) Ki(r,t) = i e-^ 1 -^ K l _ 1 {x,T)dr (* = 1,2) 

Jo 

which explicitly are given by [7]: 



~-tss-*v-0{*-v) / nf—jy- 1 



(4.3) K { = / ,/— ^ Ji-i(2y/by(t-y)) dy (i = 1, 2) 

As consequence, together with #o defined by (3.7), the other two 9 functions 

oo 

(4.4) 0i(x,t) = Yl Ki(x+2nL,t) (z = 1,2) 

n=— oo 

must be considered. 

To allow a plainer reading let's set 

(4.5) Gi{x,£,t) = 9i(\x-£\,t) + 9 i (x + £, t) (» = 0,1,2) 
In this manner, owing to (3.8) one has: 

(4.6) u(x, t) = [ [G (x,t, t) u (0 - Gxfat, t) v (0 ]d£ + 

Jo 

-2e / 6 Q (x, t-r) i/ji(t) dr + 2 e / 9o(L-x,t-r) tj} 2 (r) dr 
Jo Jo 

+ / dr I G (x,C,t-r) <p[&t, u(£,t)] d£. 
Jo Jo 

As for the v component, by (4.1) one deduces: 

(4.7) v(x,t) = voe-P* +b [ [Gi (*,£,*) u (£) -G 2 (x,£,t) Vo(t)]d£ + 



Jo 

2b e / 9\ {x, t - r) ^i(r) dr + 2 b e 9 1 (L-x,t-r) tp 2 (r) dr 
Jo Jo 

+ b I dr I Gi(x,S,t-r) p[£,r, u(£,r)]d£. 

JO JO 



Let us observe that the kernels Ki(x,t) and K<i{x,i) have the same properties 
of K (x,t). In fact [7]: 

Teorema 4.1. For all the positive constants a, be, (3 it results: 

(4.8) f \Kx\ d£ < E(t); f dr f \K t \ d£ < A 

J?R JO J?R 

(4.9) [ \K 2 (x - £,t)\ d£ < i e~ ay - p{t - y) (t-y)dy <tE(t) 
Jn Jo 

where E(t) is defined in (3.4)i and f3\ = (a/3) -1 . I 



Now, let || z || = sup | z(x, t)\, and let Bt denote the Banach space 
(4.10) B T = {z(x,t) : z £C(n T ), \\z\\ < oo }. 



By means of standard methods related to integral equations and owing to basic 
properties of K{, d (i = 0,1,2) and (p(u), it is easy to prove that the mapping 
defined by (4.6) ia a contraction of Bt in Bt and so it admits an unique fixed point 
u(x,t) G Bt [4,5]. Hence 

Teorema 4.2. When the initial data (uo,vq) are continuous functions, then the 
Neumann problem related to the non linear (FHN) system (1.2), (1.3) has a unique 
solution in the space of solutions which are regular in £It- I 



Continuous dependence for the solution of (-P/v) is an obvious consequence of 
the previous estimates. As an example of asymptotic properties let us consider the 
case tfti = ijj2 = and let 

|| (p || = sup | ip (x, t, u)\, 
n T 

then by means of (4.6), (4.7) and owing to the estimates (4.10), (3.6), (4.8), (4.9), 
the following theorem can be stated: 



Teorema 4.3. For regular solution (u,v) of the (FHN) model, whenijji = ip2 = 0, 
the following estimates hold: 



\u\ < 2 [ ||«o|| (1 + nVbt) e -wt + ||v || £(i) + (3 \\tp\ 
(4.11) 

M < \\vo\\ e-P* + 2[6(||« || + t \\v \\)E(t) + bfa. \\<p\\] 



Therefore, when t is large, the effect due to the initial disturbances (uq,vq) is 
exponentially vanishing while the effect of the non linear source is bounded for all t. 

All the previous results can be applied to the boundary Dirichelet or mixed 
conditions, too. 
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